Stability of Ge-related point defects and complexes in Ge-doped SiC>2 
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We analyze Ge-related defects in Ge-doped Si02 using first-principles density functional tech- 
niques. Ge is incorporated at the level of ~1 mol % and above. The growth conditions of Ge:Si02 
naturally set up oxygen deficiency, with vacancy concentration increasing by a factor 10 5 over un- 
doped Si02, and O vacancies binding strongly to Ge impurities. All the centers considered exhibit 
potentially EPR-active states, candidates for the identification of the Ge(n) centers. Substitutional 
Ge produces an apparent gap shrinking via its extrinsic levels. 

PACS: 61.72.Bb, 71.55.Ht, 61.72.Ji 



Silicon dioxide, besides its role in silicon-based micro- 
electronics, is the central material of fiber optics technol- 
ogy. Germanium doping of silicon dioxide is a key to a 
number of technologically relevant applications, as e.g. 
low-loss and modulated-refraction fibers. Ge can be in- 
troduced in Si02 in concentrations in the 1-10 % range 
0. As is generally true of defects in solids, the identifi- 
cation and characterization of Ge-related defect centers 
in Si02 is a difficult task, and far from completion de- 
spite intense recent investigation. This is the case, for 
instance, for the understanding of oxygen-deficient Ge- 
related centers and the so called Ge(n) [n=0, 1,2,3] cen- 
ters, for which rather disparate models have been pro- 
posed. First-principles defects theory can play a 
key role in this context, predicting the stability regimes 
and concentrations of Ge and related defects in Si02, 
and their extrinsic electronic levels and potential mag- 
netically active states. Recently, for example, these tech- 
niques were used to pinpoint the correlation of the E' 
center in Si02 with the singly-charged O vacancy 0] and 
the role of hydrogen in determining the leakage current 
across thin silica layers |J. 

Using a similar methodology, here we study the en- 
ergetics, extrinsic levels, and solubility of Ge in SiO^- 
We predict experimental signatures of a selection of 
Ge-related defects, including Ge-oxygen vacancy com- 
plexes in various stable and metastable configurations 
and charge states. We demonstrate that oxygen off- 
stoichiometry occurs naturally in Ge-doped Si02, and 
specifically that oxygen deficiency (i.e. O vacancy forma- 
tion) occurs preferentially near substitutional Gesi sites. 
The electronic structure of Gesi may explain the observed 
effective reduction of the gap in Ge-doped silica || . Sin- 
gle and paired substitutional Ge and their complexes 
with an O vacancy are all found to have accessible param- 
agnetic states: Gesi, and the GesrGesi or Gesi-Vo-Gesi 
complexes, are candidates, respectively, for the Ge(l) and 
Ge(2) electron-capturing centers; two metastable E'-like 
Ge-related centers correlate with the Ge(3) hole center. 

Our method can be summarized as follows. 
The equilibrium concentration of a defect D= 



N s exp(— Fform/keTg) is determined by the growth tem- 
perature T g , the number N s of available sites, and the 
formation free energy Ff orm =Ef orm -T Sf orm . The latter 
depends |(| on the chemical potentials of atoms added 
or removed, on the defect charge state, i.e. the charge 
released to or captured from the electrons reservoir made 
up by the whole crystal. Given the formation energies of 
the relevant defects, the concentrations and the electron 
chemical potential \i e are determined self-consistently to 
satisfy charge neutrality, as detailed in . A specific de- 
fect configuration or charge state is predicted to exist if 
its formation energy is lower than that of all other defect 
states for some value of fj, e , or when a sufficiently high 
energy barrier prevents its disappearance. The formation 
energy for a defect in charge state Q can be written as 



^form — ^dcf 



Q[i e + M(Q) - n Si ^ Si - n Gc ^ Gc - n Q MO, 



where E^ of is the total energy of the defected system 
in charge state Q, /j, e the electron chemical potential 
(equaling the Fermi level Ep in our T=0 calculations), 
and M(Q) a defect- and charge state-dependent multi- 
pole correction The nsi, iig c and no atoms con- 
sidered in the modeling of a specified defect, possess the 
chemical potentials /isi, hgc and /io, discussed further 
below. 

Energies and forces are calculated from first-principles 
within density-functional theory in the local approxi- 
mation, using the ultrasoft pseudopotential plane-wave 
method as implemented in the VASP code ||. An iso- 
lated defect is simulated in periodic boundary condi- 
tions via repeated tetragonal 71- and 72-atom super- 
cells of crystalline a-quartz Si02 having the (theoretical) 
linear dimensions 18.49, 16.02, and 20.44 atomic units. 
Atomic geometries are optimized for all Qs (which are 
modeled by removing or adding electrons as appropri- 
ate, the added charge being compensated by a uniform 
background) until all residual force components in the 
system are below 0.01 eV/A. No symmetry restriction 
is imposed. A (222) Monkhorst-Pack mesh is used || 
for k-space summation (4 points in the supercell Bril- 
louin zone). Total energy differences of different charge 
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states define ionization energies, i.e. the energy needed 
to promote e.g. an electron from the valence band into 
an empty acceptor level. Formation entropies are beyond 
the scope of the methods used here; plausible estimates 
are used when needed. 
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FIG. 1. Solubility region of Si-substituting Ge in Si02 as 
a function of the impurity and oxygen chemical potentials. 
Relevant limiting values of the latter are indicated, numbered 
as in the text. The oblique line is the solubility limit due 
to Ge02 formation. The limit given by GeO formation is 
irrelevant and is not displayed. 

To determine the formation energy of Ge and related 
defects, we must discuss the thermodynamical growth 
conditions of Ge-doped Si02- The chemical potentials 
fall within a range determined by the energy /j, hulk of the 
condensed phase of each atom ( "bulk" oxygen is assumed 
to be the oxygen molecule) , and the calculated formation 
enthalpies AH Si ° 2 =-11.31 eV of Si0 2 and AH Gc ° 2 =- 
8.80 eV of Ge02, the most stable compounds liable to 
form out of O, Si, and Ge. The following conditions ap- 
ply: 



bulk I ^jjSiC>2 



Msi 

■ bulk + 0.5AH Si ° 2 



^ ,i ..bulk 

< Msi < Msi > 
S Mo S Mo ! 



Mo 

^Si + 2^o = A*Si0 2I MGc + 2^o < MGe0 2 

The relevant extremal conditions are then 



, . , .bulk , , 

Mo — Mo > MSi — 



1. O-rich, Si- and Ge-lean 
l4f k + AH si ° 2 and ^ Gc = Mof + AH Gc ° 2 ; 

2. Ge- and Si-rich, O-lean : ^ Gc = /ij^f, m = ^f k , 



lSi0 2 



and no = Mo + °- 5 AH 

3. Ge-rich, O-lean with respect to Ge02, and Si in- 
termediate : fiGc = MGe k : Msi = Ms" lk + (AH Si ° 2 - 
AH Ge ° 2 ), and ^o = M Q ulk + 0.5 AH Gc ° 2 . 



These three cases are indicated in Fig. |[ which depicts 
the solubility region of Ge in Si02 in the {/zo,MGe} plane. 

To substitute as much Si as possible with Ge, one 
should work in comparatively Ge-rich and Si-lean con- 
ditions. Case 3 above, the most favorable in this respect, 
demands at the same time that O be lean with respect 
to Ge02 formation. As a consequence, the material will 
contain both a high Ge concentration, and a large amount 
of oxygen deficient centers (as shown below, these tend 
to localize near Ge atoms) . The minimum formation en- 
ergy of Gesi, 0.88 eV, is obtained indeed in case 3 above. 
The ensuing solubility of Ge at a tipical growth temper- 
ature T g = 1700 K (which we assume throughout) is 0.2 
mol %. This encouragingly high value falls well within 
the experimental range assuming a plausible formation 
entropy of 3-4 ke . 
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FIG. 2. Formation energy of Ge-related defects: Gesi, 
Ge-adjacent oxygen vacancy, substitutional Ge pair with and 
without O vacancy. The slope of the formation energies is 
equal to the charge state (possible values are 0, -1, and -2 in 
the present case). The vertical dotted line indicates the cal- 
culated Fermi level. All the centers possess accessible param- 
agnetic (Q—-1) states. 

Fig.|^ shows the formation energies in various charge 
states, as a function of the Fermi level and for the growth 
conditions discussed above, of all the centers considered 
here, and discussed in turn below: substitutional Gesi, 
the oxygen vacancy Vo adjacent to Gesi, the Gesi-Gesi 
pair (substituting the central Si of two neighboring tetra- 
hedra), the Gesi-Vo-Gesi complex (the oxygen bridging 
the Ge-Ge pair is missing). The dotted vertical line 
marks the Fermi level as calculated in the presence of 
all the various defects, namely Ep = E ( , + 4.65 eV at 
room temperature. We assumed a gap of 9 eV, and the 
band edge effective masses m v ~10 m e and m c ~0.3 m e 
as estimated in Ref. |J. Indeed, all the defects considered 
are in their neutral state at room temperature. 
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Gesi is neutral over most of the Fermi level range, and 
captures electrons in strong n— type conditions. Struc- 
turally, neutral Gesi is relatively trivial, with Ge-0 bonds 
(unsurprisingly slightly longer, 1.67 and 1.68 A, than the 
bulk Si-0 bonds of 1.58 and 1.59 A) and the other struc- 
tural parameters (Ge-Si distance : 3.09 A, O-Ge-0 an- 
gle : 109°, and Ge-O-Si angle : 142°) being character- 
istic of an isotropically expanded but otherwise regular 
tetrahedron. The charge state Q =— 1, one of the models 
proposed B for the paramagnetic Ge(l) center, is not 
stable in the as-grown material, but it may be observed 
in n-type or electron-irradiated material. Upon electron 
capture, the Ge substitutional moves off-center causing 
an orthorombic distortion with two short and two long 
Ge-0 bonds (differing by ~ 0.2 A), in good agreement 
with the report of Ref. @ . 

The properties of Gesi also match the main Ge-related 
optical signature in highly doped samples, namely the 
apparent reduction of the absorption gap to about 7 eV 
f|. Indeed, our calculated concentration and electronic 
structure predict an absorption into the Gesi extrinsic 
level starting around 6.5-6.8 eV with an effective final 
density of states (DOS) of about ~ 10 20 cm -3 . Assum- 
ing an impurity band width of 1 eV, likely to occur at 
such high Ge concentrations, it appears that the valence- 
to-Gesi extrinsic state is as relevant as the fundamental 
interband transition, even accounting for reduced oscilla- 
tor strength, since the effective conduction band DOS of 
Si02 is in fact about lxlO 19 cm~ 3 at room temperature. 
In summary, the large absorption red shift in Ge:Si02 is 
impurity-related, as opposed to a standard alloying ef- 
fect; a simple Sii_ 2; Ge a ;02 alloying picture would predict 
(the gap of Ge02 being 5.6 eV) a shift of less than 0.1 eV 
instead of the observed ~ 2 eV at the typical Ge concen- 
tration x ~ 0.02. The behavior just discussed is similar 
to the impurity-level (or -band) effects observed e.g. in 
InGaAsN at low N concentration |lo|) . 

We now turn to Vo-Ge centers. The isolated oxygen 
vacancy at the chemical potentials giving maximum Ge 
incorporation, has a formation energy of 3.31 eV, i.e. 0.4 
eV less than in pure Si02 in stoichiometric conditions. 
At T ff =1700 K, this implies a concentration increase of 
a mere factor of 10. A vacancy in the same simulation 
cell, but as far as possible from Gesi (about 7 A away) 
has the same formation energy, i.e. it is effectively de- 
coupled from the Ge substitution. However, when the 
vacancy and Gesi are first neighbors, the formation en- 
ergy drops to 2.0 eV (Fig. ||), implying a large attraction 
for vacancies towards Gesi centers. Now this reduced 
formation energy implies that the concentrations of Ge- 
coupled vacancies is about 10 4 that of distant vacancies, 
and 10 5 that of those in pure Si02, in quantitative agree- 
ment with the observed ratio [[ll| . Therefore the growth 
conditions of Ge-doped Si02 naturally produce an off- 
stoichiometry of oxygen: Ge-doped Si02 contains of the 
order of 10 5 more "simple" oxygen deficient centers than 



pure Si02 does. In addition this oxygen deficiency is Ge- 
biased: almost all of the vacancies are localized near Gesi 
sites, due to the attractive Gesi-vacancy effective inter- 
action, and to a lesser extent to the different solubility 
limits provided by Ge02 and Si02- This applies also, as 
discussed further below, to the vacancy bridging a pair 
of Si-substituting Ge. 

The ground state of the Ge-neighboring vacancy Vo~ 
Gesi is neutral and unpuckered - i.e. the cations neigh- 
boring the vacancy remain near the vacant site, and do 
not undergo distortions of the type involved in the E' 
center [Q. It captures electrons only for extreme n con- 
ditions. The unpuckered +1 charge state is never thermo- 
dynamically stable, which rules out, e.g., the possibility 
of associating this defect to the Ge(0) or Ge(l) centers 
|fl^2|. Instead, similarly to the undoped case Q], Vj 
becomes a metastable ground state in a puckered config- 
uration (lower than the unpuckered one by 0.15 eV) if the 
Fermi level is below about 1.6 eV ]l3[|. T hereby, it may 
be a candidate for the Ge(3) center ]lj,y,|||. The E'-like 
off-site puckering involves the Si atom adjacent to the 
vacancy, since the puckering of Ge (on the correct side of 
the vacancy Q) costs 1.3 eV more than that of Si. Conse- 
quently, the Gesi is left with an EPR-detectable unpaired 
electron, in agreement with the fact that the EPR lines of 
Ge(3) show [[l6| the characteristic signatures of hyperfine 
interaction with 73 Ge. 

Given the larger amount of vacancies and their pref- 
erence for the vicinities of Ge, and accounting for the 
73 Ge and 29 Si isotopic abundances of ~8% and 4%, the 
concentration of Si-Ge E' measured by EPR should be 
a factor of 10 4 that of E' in Si02- Observed values |l7]] 
are in a range upwards of 10 2 . In agreement with the ob- 
served axial simmetry of Ge(3), the bond lengths between 
the EPR-active Ge and the three first-neighbor oxygens 
are about the same. Also, the fact that the puckered 
configuration is lower in energy than the unpuckered one 
agrees with the observation of Ge(3) defects even in non 
irradiated samples Q. 

As mentioned, the thermodynamical stability of Si- 
Ge E' is possible only at rather extreme p conditions 
(p e ~1.6 eV), not realized in as-grown material. This 
center should therefore be observed only upon (radiation- 
induced or electrical) hole injection. In as-grown mate- 
rial, it may still be possible that, after its excitation to 
the positive state by e.g. an optical excitation, the return 
of Ge-Si E' to the neutral (and thence to the unpuckered) 
state is slowed down due to selection rules and/or disor- 
der, neither of which have been considered here. 

So far we identified possible candidates for the Ge(l) 
electron-capturing and Ge(3) E'-like centers. We now 
move on to Ge-pair defects, consisting of two corner- 
sharing Ge-centered tetrahedra. In the neutral state, the 
tetrahedra exhibit Ge-0 bonds and angles close to those 
of the isolated Gesi in their neutral state. 

As seen in Fig. 0, the formation energy of two Ge 
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impurities at neighboring substitutional sites is almost 
exactly twice that of isolated Gesi! therefore, there is 
neither a driving force for, nor an energetic hindrance 
against, the clustering of Gegi. Kinetics, and hence ther- 
mal and growth history will play the deciding role. Once 
more, the Q=-l charge state is not a ground state of the 
defect in as-grown material, but it may be observed in 
?i-doped irradiated material, and thereby become a can- 
didate for the Ge(2) paramagnetic center, as suggested 
in Rcf. Q. 

Finally we consider the formation of Vo between the 
two Ge atoms: it costs 0.15 eV less than between a Ge 
and a Si (Fig. ||), because of the previously discussed 
natural tendency of the two Ge to host off-stoichiometry 
of oxygen nearby them; this gives a concentration of a 
further factor of 3 higher than for the Ge-neighboring 
vacancy. In this case, the negative charging level is just 
slightly (~ 0.1 eV) above the Fermi level so that, given 
the uncertainties on the energy levels, it may well turn 
out to be already occupied in weakly, unintentionally 
n— type as-grown material, or will readily become occu- 
pied in moderately n— doped or irradiated material. Thus 
Ge-Vo-Ge stands out as well as a candidate for the Ge(2) 
paramagnetic center ffl. 

Once more, the unpuckered +1 vacancy is not a sta- 
ble ground state; its E'— like puckered-Ge configura- 
tion becomes metastable below /i e =1.5 eV As the 
metastable energy minimum is found to be 0.45 eV higher 
than the undistorted one, the EPR-active Ge-pair E'- 
like configuration is far less frequent than the Si-Ge one 
(~10 -7 at room temperature, although it is possible that, 
as for the standard E', a confinement barrier exists pro- 
longing its existence). This suggests that the proposed 
PI attribution of the Ge(3) center to this complex is 
unlikely to be correct. 

In summary, we have shown that thermodynamical 
growth conditions of Ge-doped Si02 naturally produce 
off-stoichiometry of oxygen, and that oxygen deficient 
centers form preferentially near Ge impurities. We find 
that energetics is neutral as to clustering of substitutional 
Ge, so that kinetics will be important. The calculated 
concentration (in agreement with experiment) and elec- 
tronic structure of Gesi show a large density of states 
starting at 6.5-6.8 eV which may explain the apparent 
gap reduction to about 7 eV in absorption. Our results 
contribute to the identification of the Ge(l), Ge(2), and 
Ge(3) centers. Ge(l) and Ge(2) are electron-capturing 
paramagnetic centers; the former is probably related to 
the singly negative states of Gesi; the latter may be as- 
sociated with the singly-negative Gesi-Gesi or Gesi-Vo- 
Gesi complexes. Ge(3) may instead be attributed to 
the Ge-Si E'-like center. Hyperfine-parameter calcula- 



tions are underway to pinpoint these attributions. In 
any case, these centers should be preferentially observed 
in n-doped or irradiated material. 
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